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Abstract. The penetrable-sphere model has been introduced in the literature to describe the peculiar thermodynamic behavior 
of some colloidal systems. In this model the interaction potential is = £ > if the two spheres are overlapped {r < a) 
and (j){r) =0 otherwise (r > C7). In this paper the shear viscosity, thermal conductivity, and self-diffusion coefficients of a 
dilute gas of penetrable spheres are evaluated. It is found that the effective collision frequency v(r*) grows as \/T* up to 
T* = ksT/e ~ 0.25, reaches a maximum at T* ~ 0.415 and then decays as r*^"'''^logr* for large temperatures. The results 
are applied to the hydrodynamic profiles in the steady Fourier and Couette flows. 

1. INTRODUCTION 

Traditionally, the kinetic theory of gases has been applied to gases made of particles which interact according to 
unbounded spherically symmetric pair potentials, such as hard spheres, power-law repulsive interactions, the square- 
well model, or the Lennard-Jones potential Moreover, granular gases are usually modeled as constituted by 
inelastic hard spheres. In the last decade, however, the equilibrium properties of fluids of particles interacting via 
bounded pair potentials have been the subject of increasing interest, the Gaussian core model and the penetrable- 
sphere model being among the most popular ones |2, 3]. These models have been proposed in the literature in 
order to understand the peculiar behavior of some colloidal systems, such as micelles in a solvent or star copolymer 
suspensions. The particles in these colloids are constituted by a small core surrounded by several attached polymeric 
arms. As a consequence of their structure, two or more of these particles allow a considerable degree of overlapping 
with a small energy cost |,2J. These are a few particular cases of systems defining what is commonly known as "soft 
matter," which has become an active field of research with interesting physical, chemical, and engineering applications. 

While the equilibrium properties and phase diagrams of particles interacting through soft potentials have been 
studied with great detail, the nonequilibrium transport properties of those systems seem to have received less attention. 
The aim of fliis paper is to contribute to the understanding of the nonequilibrium properties of particles interacting via 
bounded potentials by analyzing the transport coefficients of the penetrable-sphere (PS) model in the low-density 
regime. 

2. COLLISION PROCESS IN THE PENETRABLE- SPHERE MODEL 

The PS interaction potential is defined as 

Therefore, in this model the gas behaves as a hard-sphere (HS) gas in the low-temperature limit {T* = kgT /e 0) 
and as a collisionless gas in the high-temperature limit (T* oo). The Liouville operator and the Boltzmann-Lorentz 
collision operator for this model have been derived in Refs. |4] and |5], respectively. The PS model should not 
be confused with the Widom-Rowlinson (WR) model of interpenetrating spheres [6], which is thermodynamically 
equivalent to a binary mixture of highly non-additive hard spheres. The transport properties of the WR model are 
obtained in Ref. \J^. 

Before delving into the transport coefficients of the PS model, let us analyze its peculiar binary collision process. 
As usual, we consider the equivalent one-body problem in which a projectile particle (with the reduced mass /i =m/2 



FIGURE 1. (a) Sketch of the coUision process when ^* > 1 and b* < n{g*). The incoming incidence and refraction angles are 
a = sin^' b* and /3 — sin^' [b* /n(g*)], respectively, while the outgoing angles are /3 and a, respectively. Thus, the scattering angle 
is X = 2(/3 — a). This sketch actually corresponds to ^* = 1.2 and b* = 0.3. (b) Different possible trajectories for g* = 1.1. 



of the colHding pair) feels the central potential (j){r) centered at the origin. The projectile approaches the "target" 
with a (relative) speed g and an impact parameter b, being deflected after interaction with a scattering angle x{b* ^S*) 
that depends on the reduced impact parameter b* = b/a and reduced speed g* = g/ yjlejji. Obviously, \f g* < 1, the 
projectile does not have enough kinetic energy to "penetrate" into the core of the target and consequently it is deflected 
exactly as if the target were a hard sphere; the associated scattering angle is therefore x{b* ,g*) — 7:~2 sin^ ' b* . On the 
other hand, if g* > 1, the projectile traverses the core of the target moving with a (reduced) kinetic energy g*^ — 1 and 
eventually leaves the core along the scattering direction with the same speed g* as before collision. This penetration 
process is analogous to the double refraction of light through a sphere made of a transparent material of relative 
refraction index n{g*) = a/ 1 — g*^^ < 1. Figure^a) shows a sketch of a typical collision with ^* > 1. The scattering 
angle in this case is x{b* ^g*) = 2sin^' [/7*/n(g*)] — 2sin^' /?*. However, the latter expression is valid only if the impact 
parameter is sufficiently small, namely iib* <n {g* ) . Otherwise, the incidence angle a = sin^ ' /?* is large enough as to 
produce a "total reflection" effect, so the projectile is again deflected as if colliding with a hard sphere [cf. Fig.^b)]. 
In summary, the scattering angle xi^* tS*) for the PS potential is given by 



cosx{b*,g*)=2b*^ -\ 



Mb*:8*), g* >landO<b* <n{g*), 
0, otherwise. 



(2) 



where 



Ai{b*,g*)^2Jl-b*^/n 



. (b*^/n{g*)^ ^\-b*^- {ib*^ - l) ^\-b*^/n^(g*) 



(3) 



At a given value of the (reduced) relative speed > 1, Z starts increasing with b* until it reaches a maximum value 
Zmax(§*) = cos^' (l — 2g*^^) at b* — n{g*), and decreases thereafter. Figure|2]shows co?,x for a few cases. 

From the scattering law one can obtain the differential cross section as B{x,g*) = j C7^| (9/7*^/(9 cos If g* < 1, 
then B{x,g*) = ic^, as for hard spheres. Otherwise, 



Bix,8*)/cy' = 



0, 



b*<n{g*) 



X > Zmaxfe*), 
, 0<X<Xm!,xU 



(4) 



where 

db*^ 



dcosx 



2Vl+cosx(2^*2 _ 1)+^* ^2(^*^-1) [l +3C0SX -2(1 - cosx f 8*^8*^ - 1) 



b*<n{g*) 



4^\+cosx[l+2{l-cosx)8*^{g*^ - 1)]^ 



(5) 

Figure |3] shows the differential cross section B{x,g*) for several values of ^* > 1. We observe that B{x,g*) =0 for 



X > Zmaxte*), i.e. for cos;i; <l-2g* \ takes the HS value B{x,g*) = |c7^ x X, 



and monotonically 




increases with decreasing scattering angle for % < Xmaxig*), i-C- for cosj > 1 — 2g*^^. This imphes that, as expected, 
when the kinetic energy is much larger than £ most of the collisions are grazing. On the other hand, the total cross 
section, which is given by the area below the curves in Fig.|3j is independent of g* and hence it coincides with that of 
HS, namely In Jl_y ,2d{cosx)B{x,8*) ^ na^. 



3. TRANSPORT COEFFICIENTS 

As is well known, the Chapman-Enskog method allows one to derive the Navier-Stokes transport coefficients from 
the Boltzmann equation for a dilute gas in terms of the scattering law corresponding to the interaction potential of 
interest 1 1]. For a single gas, the relevant transport coefficients are the shear viscosity Tj(r), the thermal conductivity 
k{T), and the self-diffusion coefficient D{T). In the first Sonine approximation, their expressions are QJ 

, , 5 kijT , , 15 fen , , , , 3 kaT 

^^^877^' ^iT)^^-ri{T), D{T)^- ^ (6) 

where the coUisional integrals £lk/{T) are 



^^,,(7-) = ^^f^'y^'^'y^'^'Qi {y^2k^) , Qt{g) = ^^[dbb [l -cos^;t( 
In the special case of hard spheres, one has 



(7) 



r,HS,7,^ 2 (^+1)' ^ 2(1+^) 



5 y^mkBT/n 15 feg 3 y^ksT/mK 

16 4 m 8 n(7^ 

Now we particularize to penetrable spheres. In that case, it is convenient to define the reduced integrals Q.l^{T* ) = 
Q.k,e{T)/Q.ff{T), so that 

77hs/j?= 0^2(7'*), Kns/K^al2iT*), Dns/D^al,iT*). (10) 
The functions Q.lf{T*) can be expressed as 

^lAT*) = 1 r ^ ,y 1 r dye-y'y^'+'R, (yVf^) , (11) 



2(1+^) 



(A:+l)! 
where we have called 

Rdg*)^l"^' ^ db*b*Aeib*,g*), Ai{b*,g*)=cos^xib*,g*)-{2b*^~iy . (12) 
In particular, taking into account Eq. and setting £ ~ I and £ ~2, one gets 



g*^-l)(g*^+2)+-^(g*^-iy''\ 4g*- - 4g* 



,*4 



*2-l) (2g*^-l)ln(2g*^~2g*Jg*^-l-l\, (13) 



1 / .9 .\-2 



-,*6 



2g*^-l)(2g*'~2g*^ + l]ln[2g*'~2g*Jg*'-l-l]. (14) 



Insertion of Eqs. (I13> and J14> into Eq. Jl 1> gives i2[ j (T*) and ilj 2(7^*) in terms of quadratures that need to be 
evaluated numerically. On the other hand, for high temperatures (T* » 1) the integral in Eq. (II 11 1 is dominated by the 
behavior of Re{g*) for g* ^ 1. In that limit Eqs. jl3l l and J14l i become 

„,„.)4-i±|iP^)+o(«-,„s-), w)4-i±|ie^)+o(«-«,„s-)^ ,15, 

This yields the following asymptotic approximations: 

i2i,(r)«i-.-vr^i + i_)_!i(_±p, (16) 
2(r ) « 1 - f 1 + ± - 'J^) - ^'(-^f), (17) 

where Ei(z) = — r(0, — z) = — J^.dtt^^e^' is the exponential integral function, r(n,z) being the incomplete gamma 
function. Quite interestingly, the approximate expressions J16t and illl agree almost perfectly for the whole temper- 
ature domain with the results obtained from Eq. il l\ by numerical integration, as shown in the left panel of Fig.|3 As 
usual, we can introduce an effective collision frequency associated with the shear viscosity as (ill v(r) = nksT /rjiT). 
Therefore, 

v(r*) = vo\/7^i2^.2(r*), vo = ynaV^e/'w- (18) 

In the case of hard spheres, Vus{T*) = VqVt*- The right panel of Fig.lHshows v(T*)/vo for PS and US. We observe 
that up to T* ~ 0.25 both systems have practically the same collision frequency. For larger temperatures, however, both 
collision frequencies strongly differ each other. While in the HS case the mean free path is independent of temperature 
and hence Vhs(7^*) grows proportionally to the thermal velocity, penetrability effects in the PS model become more 
and more important for T* > 0.25. As a consequence, v{T*) reaches a maximum value Vmax — 0.548vo at T* = T^^x — 

0.415, and then decays for asymptotically large temperatures as v{T*) w ^^jjiT* + j — ^7+ln-\/2^, where 

7~ 0.577216 is Euler's constant. 
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FIGURE 4. Left panel: £1* [(r*) and Q.22{T*) as obtained from Eq. <1 U by numerical integration (solid lines) and from the 
approximate expressions <16> and <17t (dashed lines). The numerical and approximate curves are practically indistinguishable in 
the case of flj 2(^*)- Right panel: Effective collision frequency v(r*)/vo for the PS model (solid line) and for HS (dashed line). 



4. APPLICATION TO THE FOURIER AND COUETTE FLOWS 

Let us compare now the temperature and velocity profiles for HS and PS in the steady Fourier and Couette flows 
10]. In the planar Fourier flow the gas is enclosed between two parallel plates at rest located at y — and y — L 
and kept at temperatures T\ and 72, respectively. In the steady state described by the Navier-Stokes equations the 
temperature profile is the solution to v^^ {T)dT /dy = const. Neglecting boundary layer effects and applying the 
boundary conditions, the implicit solution is 



y/L. 



T* 



dev-\e) 



-1 



r* 



dev-\e). 



(19) 



In the case of hard spheres, this yields T* (y) = T* [l + {y^TfJlf ~ l)y/L\ . 

In the planar Couette flow the plates aty — Q and y — L move along the jc-direction with velocities —U and +U, 
respectively, but otherwise they are kept at the same temperature Tq. According to the Navier-Stokes equation in the 



To* [1 +m2(1 . 



where 



steady state the temperature and velocity fields are related |8, pp. 213-222] by T*{u* 
u* —Ux/U is the reduced flow velocity and M = ^JlmU^ /ISksTQ is a sort of Mach number. The velocity profile obeys 
the condition v^^ [T {ux))dux / dy — const, whose implicit solution (again neglecting boundary effects) is 



j^^dwv-\T*{w)) p dwv-\T*{w)). 



(20) 



For hard spheres, ul(y) = Vl+M-^sin [(2y/L- l)sin"'(l +M-2)-i/2] 
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FIGURE 5. Left panel: Temperature profiles for (a) the Fourier flow with Tj* =0.1 and T2 = 1 and (b) the Couette flow with 
To = 0. 1 and M = ^/imU^JTSk^) = 3. Right panel: Velocity profile for the Couette flow with Tq = 0. 1 and M = 3. 



The temperature and velocity profiles corresponding to HS and PS are compared in Fig. |5] for representative 
examples of Fourier and Couette flows. Note that in the region where T* > 7^^^ — 0.415, the PS and HS systems 
have opposite curvatures for the temperature (Fourier) and velocity (Couette) profiles. 



5. CONCLUDING REMARKS 

The aim of this paper has been to derive and analyze the Navier-Stokes transport coefficients of a dilute gas made 
of penetrable spheres (PS). Apart from its ability to describe qualitatively the effective interactions in some colloidal 
systems, the PS model is also interesting as a crossover between a gas of hard spheres (HS) in the low temperature 
limit (r* 0) and a gas of non-interacting particles in the opposite limit {T* °o). Although the total cross section 
of PS at any temperature is the same as that of HS, the effective collision frequency v{T*) associated with the shear 
viscosity and thermal conductivity exhibits a non-trivial temperature dependence. As expected on physical grounds, 
v(T) ranges from v{T*) VusiT*) °^ Vf* when T* ^ to v{T*) when T* oo, so it presents a maximum 
value Vmax at a certain temperature Ij^ax- On the other hand, plausible physical arguments do not suffice to answer 
more specific questions, for instance: (i) at which threshold temperature T^Jj does v{T*) start to deviate significantly 
from Vhs (ii) what is the value of T^^^J; (iii) how small is Vmax as compared with VHs(7"max)^' (i^) ^'^ which 

temperature T^^2 > ^max the collision frequency has decayed to half its maximum value?; (v) how is the asymptotic 
behavior of v(r*) in the high-temperature limit? 

The evaluation of the transport coefficients and of v{T*) has required the detailed analysis of the collision process. 
In principle, the collision integrals must be evaluated numerically [cf. Eqs. il 1I I-J14H. but excellent approximations 
are provided by Eqs. JTEi and ( [T7t . The answer to the questions posed above are: (i) T^^ ~ 0.25; (ii) T^^^ ~ 0.415; 
(iii) Vmax/vHs(7'max) - 0-850; (iv) r*/2 ^ 1.619; (v) v{T*) « ir*"^/^lnr* for T* oo. The conti-asting temperature 
dependencies of the PS and HS transport coefficients for T* > 7[[^ ~ 0.25 have been illustrated here by comparing the 
respective hydrodynamic profiles for the Fourier and Couette flows. 
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